Electron Dynamics in a DNA Molecule 
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We report on our theoretical investigations of the electronic states in a DNA molecule. We have 
used a two-leg charge ladder model where electron-electron interactions and the electron spin have 
been taken into account. The energy spectra for G-C and A-T base pairs obtained by numerically 
diagonalizing the Hamiltonian reveal a gap structure and the interaction is found to enhance the 
energy gaps. We also present the charge distribution in the ground state and low-lying excited states 
for the A-T and G-C base pairs. 



The unique properties of DNA, self-assembly and 
molecular recognition, has rendered the 'molecule of life' 
[1] a promising candidate in the rapidly emerging field 
of molecular nano-electronics [2]. A recent report of a 
field-effect transistor based on DNA molecules [3], that 
was preceded by a series of seminal experiments on the 
electron conduction in DNA [4], has sparked a lot of in- 
terest on the electronic properties of the DNA. A thor- 
ough understanding of the electronic properties of DNA 
is crucial in the development of the future DNA-based 
nanoscalc devices. In addition, charge transfer through 
DNA also plays an important role in radiation damage 
and repair [5] and therefore important for biological pro- 
cesses. Following the techniques pioneered earlier for 
a direct measurement of electrical transport through a 
DNA molecule [6], Yoo et al. [3] measured the conduc- 
tance of poly(dG)-poly(dC) DNA (1.7 - 2.9 ^m long) 
and poly(dA)-poly(dT) DNA (0.5 - 1.5 /xm). The bundle 
of DNA molecules were trapped between electrodes that 
were 20 nm apart. The I-V results showed a strong tem- 
perature dependence of the current and was interpreted 
in terms of a small polaron hopping model where the cur- 
rent was expressed as I oc sinh/3V e _Ba ' feT . Here, E a is 
the activation energy, T is the background temperature, 
[3 = ea/2kTd, e is the electron charge, a is the hopping 
distance, and d is the separation between the electrodes. 
The results for poly(dA)-poly(dT) DNA can very well 
be described in this way if f3 is taken to be indepen- 
dent of temperature. In the case of poly(dG)-poly(dC) 
DNA molecules, a similar temperature dependence was 
observed but with a much lower resistance of 1.3 Mf2 
at room temperature as compared to 100 Mil for the 
P oly(dA)-poly(dT) DNA. The P oly(dG)-poly(dC) DNA 
also shows the correct temperature dependence of (3 and 
the temperature dependence of the current was observed 
down to 4.2 K. In contrast, for the poly(dA)-poly(dT) 
DNA, temperature dependence of the current was ob- 
served down only to 50 K. 

Yoo et al. also performed I-V measurements with 
an additional electric field from the back of the Si sub- 
strate and measured the FET property of the back-gate 
type. They noted that, in the FET based on poly(dA)- 
poly(dT), application of a positive gate voltage results 



in a larger conductance, i.e., an indication of n-type be- 
havior. In the case of the DNA-FET based on poly(dG)- 
poly(dC), a p-type conducting behavior was observed. 

The DNA conductivity measurements discussed above 
have also inspired a few theoretical studies [7-14]. These 
works primarily focused on the evaluation of transport 
through one-dimensional systems using a model Hamilto- 
nian, or electronic structure calculations from first princi- 
ples. It should be pointed out that the effect of electron- 
electron interactions is important for understanding the 
physical properties of the DNA molecules, such as exci- 
tation spectra, charge distribution and charge transport 
in DNA molecules. The reason for that is the following: 
There are a few energy scales which determine the DNA 
properties. The first one is the tunneling gap, or hopping 
integrals between nearest DNA base pairs. The typical 
value of these hopping integrals is 0.5 eV. The second 
energy scale is determined by the single-particle energy 
spectrum of a single DNA base pair. This is the energy 
between the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO), 
which is of the order of 1 eV [7]. The other energy scale 
comes from interactions between electrons at different 
and at the same base pairs. This energy scale is 1 eV 
for Hartree interactions and since it is close to HOMO to 
LUMO excitation energy the interaction can have strong 
effect on many-particle excitation spectra and charge dis- 
tribution. At the same time the exchange interaction, 
which can be of the order of hopping integrals between 
nearest base pairs, should have strong effect on electron 
transport. 

In this paper, we report on the electronic properties of 
the DNA, in particular, the influence of electron-electron 
interaction on the energy spectrum and the excitation 
gap. We model the double-stranded DNA as a two-leg 
charge ladder [7]. As a first approximation, we con- 
sider only the electronic degree of freedom and disre- 
gard the vibrational modes and their effects on the elec- 
tronic motion [15]. The Hamiltonian of the electronic 
system consists of two parts: the tight-binding Hamilto- 
nian, Ti. t , and interaction Hamiltonian, Tii. The tight- 
binding Hamiltonian is a two-chain Anderson Hamilto- 
nian describing the hopping between the nearest neigh- 
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bors (nearest bases) within each chain and the hopping 
between the two chains (within each base) 

H t = £ha i* a h* + e « 6 l A,* + 
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where e/j, e; are the (site) energies of the HOMO and 
LUMO for a single isolated base pair respectively, aj i(T , 
&i j(T are annihilation operators of electron with spin a in 
HOMO and LUMO states of i-th base pair, th is the hop- 
ping integral between HOMO of the nearest base pairs, 
ti is the hopping integral between LUMO of the near- 
est base pairs, and thi is the hopping integral between 
strands (HOMO and LUMO) of the same base pair. The 
index i labels the different base pairs, while a =\\ is the 
spin index. The summation over index i goes from 1 to 
N where N is the number of base pairs. 

The interaction part of the Hamiltonian has the fol- 
lowing form 
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which is described by ten parameters. 

In the ground state the number of electrons is equal to 
IN ', so that all HOMO states (with both spin directions) 
are occupied. To find the excitation gap and the energy 
spectrum of the electron system with the Hamiltonian 
Eqs. (1,2), we consider the DNA structure to have a finite 
number TV of base pairs and, by exactly diagonalizing 
the Hamiltonian matrix, we obtain the ground state and 



TABLE I. Paramaters of the DNA structure used in our 
present work. Energies are in eV. The subscripts h and I 
correspond to the HOMO and LUMO. 



Hopping integrals 
[16,17] 


G-C 


t h = -0.1419, t t = 0.0525 


A T 


ih — — u.uuyj, Li — u.iuj^ 


The site energies [9] 


Li" V 


£h — — 14./ 14, £[ — — I0.0U0 


A-T 


e h = -14.635, ei = -13.734 


On-site 

interaction [18] 


G-C 


V M = 5.879, Vio = 5.227 


A-T 


V h0 = 5.681, V l0 = 5.23 


Interstrand 
interaction [9] 


G-C 


Vhi = 1.844, V n = 2.455 


A-T 


Vhi = 1.625, V n = 2.378 


Intrastrand 
interaction [9] 


G-C 


Vhio = 2.7 


A-T 


Vmo = 2.6 



the lowest excitation states of the system. Some of the 
parameters of the DNA structure relevant for our studies 
are listed in Table I. To eliminate the effects of boundaries 
we have also imposed the periodic boundary conditions, 
so that in Eqs. (l)-(2) we have cin+i = &i and b^v+i = &i- 

We have considered a segment of DNA containing ten 
base pairs of poly(dG)-poly(dC) and poly(dA)-poly(dT) 
DNA molecules. It should be pointed out that in the 
experiments of Porath et al. [6], the DNA sample has 
30 base pairs. Similarly, we can estimate that since the 
distance between electrodes in the experiments by Yoo et 
al. [3] was about 20 nm, taking into account the fact that 
the distance between the base pairs is 0.34 nm, there were 
~ 50 base pairs between electrodes of the set up used by 
Yoo et al. However, in our present model, handling more 
than 10 base pairs would be a formidable endeavor and 
has not been attempted. 

From the exact diagonalization of the Hamiltonian we 
have obtained the energy spectra, E n , and the corre- 
sponding wave functions, of the DNA system. The 
energy spectra for the A-T and G-C base pairs are shown 
in Fig. 1. From this figure it is clear that the excita- 
tion gap for the G-C base pairs, Ag-c ~ L12 eV, is 
larger than for the A-T base pairs, Aa-t ~ 0.79 eV. In 
both cases the lowest excitation does not include elec- 
tron spin-flip. The spin-flip excitation (open circles) gap 
is equal to Ao-Cspin ~ 1-17 eV for the G-C base pairs 
and AA-T,spin ~ 0-94 eV for the A-T pairs. The main 
difference between the A-T and G-C pairs is that for the 
G-C pairs the lowest spin-flip excitations are very close to 
the excitations without spin reversal. For example, for 
the G-C pairs the energy difference between these two 
types of excitations is Ac-c,spin — Aq-c ~ 0.05 eV, 
which is smaller than corresponding value for the A-T 
pairs, A A - T . spin - A^-t ~ 0.15 eV. 

It is interesting to compare the energy gaps of the in- 
teracting DNA system to the corresponding gaps of the 
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FIG. 1. The energy spectra of the A-T and the G-C base 
pairs in our model. The closed circles correspond to the sys- 
tem with equal number of up- and down-spin electrons. The 
open circles correspond to the spin-flip excitations. 



non-interacting system. It follows from Eq. (1) that the 
energy gaps of the non-interacting system is equal to 
A = ei - e h - 2(ti - t h ), which gives A ,g-c ~ 1-02 
cV for the G-C pairs and A ,a-t « 0.56 eV for the A-T 
pairs. Comparing these values to the energy gaps of cor- 
responding interacting system we conclude that both for 
G-C and for the A-T pairs we have interaction enhance- 
ment of the energy gaps by 0.1 eV (for the G-C pairs) 
and by 0.23 eV (for the A-T pairs). Again we see that 
the effect of interaction is more pronounced for the A-T 
pairs than for the G-C pairs. From all these results we 
conclude that the interaction has a weaker effect on the 
energy spectra for the G-C base pairs than for the A-T 
pairs. 

In Ref. [3] the conductance of DNA molecules was 
found to have activated dependence on temperature. Ac- 
tivation energies were extracted in Ref. [3] to be 0.18 eV 
for A-T base pairs and 0.12 eV for G-C pairs at high 
temperatures. These activation energies are close to our 
interaction enhancement of the energy gaps as described 
above. The activated nature of the DNA conductance 
means that electrons during their transport through a 
DNA molecule must overcome some potential barrier. 
Our results indicate that inter-electron interactions have 
considerable contribution to the activation energy of the 
electron transport. 

The numerically generated wave functions allow us to 
calculate the electron charge distribution along the DNA 
molecule from the equations 

Pn,]{k) = ^ ^ \^ n (i 1 ,...,i Nr ;j 1 ,...,j Nl )\ 2 
x 6 (k - ii) , 
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FIG. 2. The density of electrons with up- and down-spins 
is shown for the A-T base pairs for (a) the ground state; (b) 
the first excited state with equal number of up- and down-spin 
electrons; (c) the first spin-flip excited state. 



h ■••ijv T ji---]N i 

x 5(k- ji) , 



(3) 



where p n ,a(k) is the density of the electrons with spin 
a =| or I in the state n at the base pair k, and 
^ n . . . ,ijv T ;ji, • • ■ jj'iVi) is the wave function of the 
state n with electrons with spin a =| and N± elec- 
trons with spin a =J.. Here ii, . . . , ijv T and j\, . . . , are 
the coordinates of electrons in the base-pair representa- 
tion. The results for the charge density of the different 
base pairs are shown in Fig. 2 (for A-T) and in Fig. 3 
(for G-C). 

Clearly, for the A-T and the G-C pairs the charge dis- 
tribution in the ground state is the same, which indicates 
that the electrons are strongly localized at the base pairs, 
i.e. at each base pair the electron density is equal to ei- 
ther or 1 (the difference from or 1 is less than 0.01). 
However, the charge distribution in the excited states 
of G-C and A-T systems show a different behavior (see 
Fig. 2b and Fig. 3b). For the G-C pairs, distribution 
of electrons with spin a =\ (Fig. 3b) is the same as in 
the ground state (Fig. 3a). The excitation manifests it- 
self only in the redistribution of the electrons with spin 
u =|, making the single-electron states more delocalized. 
For the A-T pairs the charge distribution in the excited 
state is different from the ground state for both spin di- 
rections (Fig. 2b). Contrary to the case of the G-C pairs, 
here the single-electron states remain strongly localized 
even in the excited state. 

The spin-flip excitations for the G-C and the A-T pairs 
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FIG. 3. Same as in Fig. 2 but for the G-C base pairs. 



also behave differently. For the A-T pairs the excita- 
tion is just spin-flip at one base (with base index i = 2 
in Fig. 2c) with a weak redistribution (delocalization) of 
other electrons over the base pairs. For the G-C pairs the 
spin-flip excitation corresponds to spin-flip at one base 
pair with index i = 2 (see Fig. 3c) with subsequent hop- 
ping of electron with spin a = J. from base pair i = 3 to the 
base pair i = 2. Therefore, in the case of the G-C pairs 
the many-particle spin-flip excitation is a single-particle 
spin-flip hopping excitation. Since the electron-electron 
interactions tend to suppress the hopping processes, the 
difference in spin-flip excitations for the G-C and the A-T 
pairs also illustrate that the effect of electron interactions 
is less pronounced for the G-C base pairs than for the A- 
T pairs. 

In summary, we have performed theoretical calcula- 
tions of the electron energy spectrum, based on a two- 
leg charge ladder model for the poly(dA)-poly(dT) DNA 
and poly(dG)-poly(dC) DNA molecules. We take the 
electron-electron interactions and the electron spin de- 
gree of freedom fully into account in our model. The 
energy spectra for the G-C and the A-T base pairs show 
a large gap and the interaction was found to enhance the 
gap. The effect of interaction is less pronounced for the 
G-C base pairs than that of the A-T pairs. The spin- 
flip excitations are not the lowest energy excitations. We 
also analyze the charge distribution for the ground state 
as well as for the excitations. The present report is the 
first step in our investigation of the electronic properties 
of the DNA. In our calculations of the energy spectra we 
have not included the vibrational modes. These modes 
are very soft in the DNA and can have strong effects on 
the excitation spectra, resulting in polaronic effects and 



strong renormalization of electron-electron interactions. 
Such renormalization of inter-electron interactions should 
in turn have strong dependence on the temperature. This 
is because the temperature, inducing the excitations of 
soft vibrational modes, strongly affect the distances be- 
tween the electrons and as a result modify inter-electron 
interactions. The effect of vibrational modes on the en- 
ergy spectra and charge distribution of DNA molecules, 
as well as systems containing more base pairs will be the 
subject of our future works. 
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